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I.  INTRODUCTION 

1-3 

An  interesting  electromagnetic  radiation  mechanism  known  as 
the  electron  cyclotron  maser  has  been  the  subject  of  intense  research 
activities  in  recent  years.  The  main  motivation  behind  these  efforts 
is  to  use  this  mechanism  in  a new  class  of  microwave  devices  called 
gyrotrons  capable  of  generating  microwaves  at  unprecendented  power 
levels  at  centimeter  through  submillimeter  wavelengths.  As  in  most 
microwave  devices,  the  basic  mechanism  for  radiation  in  the  gyrotron 
is  electron  bunching  in  the  presence  of  the  wave  fields.  To  illustrate 
this  process  we  consider  two  electrons  in  an  external  magnetic  field 
B„  e and  an  oscillating  wave  electric  field  E e (=  E cos  wte  ). 

We  assume  that  the  two  electrons  have  the  same  energy,  hence  they 
rotate  on  the  same  cyclotron  orbit  with  the  frequency 

e B^/ymc  (1) 

where  ym  is  the  relativistic  mass  of  the  electrons.  Figure  1 shows 
the  instantaneous  positions  of  the  two  electrons  at  t=0  and  the 
orientations  of  B and  E . The  two  electrons  are  initially  separated 

O X 

by  180°  in  the  phase  angle.  Accelerated  by  the  electric  field, 
electron  1 will  gain  energy,  become  heavier,  and  hence  rotate  slower. 
Electron  2,  however,  will  lose  energy,  become  lighter,  and  hence 
rotate  faster.  As  a result,  the  two  electrons  will  approach  each 
other  in  their  common  cyclotron  orbit  indicating  bunching.  In  contrast 
to  usual  bunching  mechanisms,  here  an  electron  which  is  accelerated  will 
actually  slip  back  in  phase  angle  as  though  it  had  a negative  mass. 

Hence  this  bunching  phenomenon  can  be  viewed  as  a negative  mass  effect. 
To  show  how  energy  can  be  transferred  from  the  electrons  to  the  wave, 
we  picture  that  a large  number  of  electrons  are  uniformly  distributed 
along  the  cyclotron  orbit  shown  in  Figure  1.  The  rate  of  total  electron 
energy  gain  is 

Note;  Menuacript  tubmltted  June  6,  1978. 
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R = -Me  Eq  cos  u)t. 


(2) 


where  is  the  x-component  of  the  center  of  mass  velocity  of  all 
electrons.  Initially,  R=0  because  the  center  of  mass  is  stationary. 

However,  the  negative  mass  effect  will  subsequently  cause  the  electrons 
to  bunch  so  that  their  center  of  mass  will  shift  to  the  right  and 
start  to  rotate  counterclockwise  at  the  cyclotron  frequency  Thus 
we  may  approximately  write  v^^  as 

''cx  = -''o^^^  -C^’ 

' 

where  Vp(t)  is  a slowly  varying  function  of  t.  One  notes  that  even  in 

the  absence  of  the  negative  mass  effect,  the  center  of  mass  would  ; 

oscillate  in  the  presence  of  the  wave  electric  field.  An  exact  expression 

for  Vex  should  include  this  oscillation  as  additional  terms  in  Eq.  (3). 

However,  to  simplify  this  example,  we  have  only  considered  the  negative 
mass  contribution  here.  In  the  formal  analysis  which  follows,  the 
complete  electron  dynamics  are  included.  Inserting  Eq.  (3)  into 
Eq.  (2),  we  obtain 

R = Ne  VQ(t)  Eq  sin  S^^t  cos  wt 
= Ne  VpCt)  Eq  [sin(f2Q-u))t  + sin(fiQ+a))t].  (4) 


If  Qq  u),  then  the  first  and  second  terms  in  Eq.  (4)  are,  respectively, 
a slowly  and  rapidly  varying  function  of  t.  Thus,  on  a slow  time  scale 
(i.e.,  neglect  fast  oscillations),  Eq.  (4)  can  be  written 

R s Ne  VQ(t)  Eq  sin(i2Q-a))t.  (5) 
From  Eq.  (5),  we  find  that  negative  electron  energy  gain  occurs  if 
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,Vv  ' 


(6) 


Brief  summaries  of  gyrotron  theories  and  experiernents 
together  with  lists  of  references  can  be  found  in  recent  review 
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papers.  ’ Reference  4 also  contains  a good  summary  of  various  gyrotron 
configurations.  In  the  present  study,  we  will  concentrate  on  a 
particular  type  of  gyrotron  configurations  - the  gyrotron  traveling 
wave  amplifier  (gyro-TWA)  which  consists  of  an  electron  beam  propagating 
in  a simple  waveguide.  As  compared  with  the  other  types  of  gyrotrons. 
gyromonotron  and  gyroklystron,  which  employ  resonant  cavities,  the 
gyro-TWA  generally  has  broader  bandwidth  and  higher  power  handling 
capability.  The  operation  of  a gyro-TWA  was  first  demonstrated  by 
Granatstein,^  et  al.,  using  an  intense  relativistic  electron  beam.  A 
gain  of  16  dB  (1.1  dB/cm)  was  achieved  at  8.6  GHz  with  a bandwidth 
^5%  and  an  output  power  of  4 MW.  However,  efficiency  was  relatively 
low  (<1%). 

Cyclotron  maser  interactions  in  a gyro-TWA  have  generally  been 

treated  as  an  electromagnetic  instability  driven  by  the  velocity 

anisotropy  of  the  electron  medium. Linear  theories  can  be  found 

in  Refs.  9-12,  while  nonlinear  saturation  mechanisms  have  been  studied 

in  detail  in  Refs.  7 and  8.  While  the  fundamental  physical  processes 

have  been  carefully  analyzed  in  these  papers,  the  problem  of  parameter 

optimization,  which  is  of  vital  importance  to  the  actual  operation  of 

gyro-TWA' s,  has  not  been  adequately  addressed.  This  can  be  seen  from 

the  fact  that  most  theories^’^’^®’^^  have  been  formulated  in  Cartesian 

geometry  rather  than  the  realistic  cylindrical  geometry.  In  order  to 

13 

bridge  this  gap  between  theories  and  experiments,  in  a recent  paper 
we  have  reformulated  the  linear  theory  and  developed  a numerical 

Q 

simulation  code  in  cylindrical  geometry  to  study  the  operation  of 
the  gyro-TWA  at  the  fundamental  cyclotron  harmonic.  There  the  proce- 
dures of  parameter  optimization  were  detailed  and  extensive  numerical 
data  generated.  The  present  paper  contains  a similar  study  for 
operations  at  the  nonfundamental  cyclotron  harmonics.  It  presents 


for  the  first  time  a comprehensive  simulation  study  of  the  nonfundamental 
harmonic  interactions  in  a gyro-TWA,  Also  analytical  scaling  relations 
for  growth  rate  and  efficiency  have  been  derived  to  compare  with 
the  simulation  results  and  to  scale  to  parameter  regimes  not  covered  by 
the  simulations. 

In  Section  II,  we  present  the  linear  analysis.  In  Section  III,  an 
analytical  scaling  relation  for  the  efficiency  is  derived.  In  Section  IV, 
we  discuss  the  nonlinear  saturation  mechanisms  and  present  the  results 
of  single  wave  numerical  simulations.  Section  V applies  the  present 
theory  to  the  design  of  a 94  GHz  gyrotron  traveling  wave  amplifier  based  on 
nonfundamental  harmonic  interactions.  In  Section  IV,  the  limitations 
of  the  present  theory  and  the  difficulties  involved  in  actual  experiments 
are  considered. 


i.^4  Wi 


II. 


LINEAR  THEORY 


A.  Derivation  of  the  Dispersion  Relation 

The  typical  configuration  of  a gyro-TWA  consists  of  an  annular 

electron  beam  propagating  inside  a waveguide  of  circular  cross  section 

(of  radius  r ).  The  electrons,  guided  by  a uniform  magnetic  field 
w 

(B^  e^),  move  along  helical  trajectories.  The  projections  of  these 
trajectories  on  the  cross-sectional  plane  are  shown  in  Figure  2. 

The  following  simplifying  assumptions  are  made  for  both  the 
linear  theory  and  the  simulation: 

1.  The  electrons  are  initially  monoenergetic  with 
their  guiding  centers  uniformly  distributed  on 

a surface  of  constant  radius  r . 

0 

2.  The  beam  is  sufficiently  tenuous  that  its  space 
charge  electric  field  can  be  neglected  and  the 
spatial  structure  of  the  vacuum  waveguide  mode  is 
unaffected  by  the  presence  of  the  beam.^'+ 

3.  The  beam  interacts  with  the  azimuthal ly  s^metric  TE^^ 
waveguide  mode  so  that  all  variables  are  independent 
of  the  azimuthal  angle  0. 

An  important  feature  of  the  present  model  is  the  treatment  of  the 

electron  guiding  center  radius  r as  a variable  parameter  (in  earlier 

7 ft  1 f)  ^ 

parallel-plate  models,  ’’  the  guiding  centers  were  centrally  located 
between  the  plates).  This  gives  us  the  means  to  analytically  optimize 
the  beam-wave  coupling  with  respect  to  r^  (see  item  D).  the  electron 
Larmor  radius  has  also  been  kept  as  a free  parameter  with  no  approxima- 
tion. This  is  necessary  since  nonfundamental  cyclotron  harmonic 
interactions  result  from  finite  Larmor  radius  effects. 


The  dynamics  of  the  electron  beam  is  described  by  the  linearized 
form  of  the  relativistic  Vlasov  equation. 


M ^ I ' U ' c ^ ^0  ~z  ' 9^P 


:(1) 


= e 


(E<‘>  * V X B<‘>)  ■ ^ f„ 


(7) 


where  and  are  the  initial  and  perturbed  distribution  functions, 

respectively,  and  are  the  wave  fields  of  the  TE  waveguide 


mode,  i .e. , 


(1)  = 


(8) 


(1) 


B = B,  Re  [-ik  a exp  (ik  z-i(ijt)]  , (9) 

r z znin  z 


and 


(1)  _ 


B^^^^  / k^  c . 


(10) 


In  Eqs.  (8)  thru  (10),  B^  is  the  amplitude  of  the  wave  magnetic  field, 

fx  "X  / r , and  x„  is  the  n-th  root  of  J,(x)  = 0.  Note  that  the 
n n ' w’  n ^ l /jj 

validity  of  the  linear  theory  requires  B^  B^  and  f « f^.  From 

Eq.  (8)  and  the  Maxwell  equations,  we  obtain 


iJ/  c^  - k/  - a/)  B^jQ(a^r) 


= - ^ exp  (-ik,  z + iwt)  7 ^ (r 
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r 3r 


(11) 


where  is  the  perturbed  beam  current  given  by 


(1) 


=/' 


= -e  I f^^^  Vg  d^  p . 


m m pi 


\r  ^ 
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ti 


(12) 


Multiplying  Eq.  (11)  by  rj^(u^r)  and  integrating  over  r from  o to  r^, 
we  obtain 


-Siru  exp(-ik_  2+iwt) 


cr  ^(x  ) 
wo  '■'n^ 


/ 


rJg^^^J^(a^r)  dr  , 


(13) 


Equations  (7),  (12),  and  (13)  form  a set  of  self  consistent  equations. 

To  derive  the  dispersion  relation  from  these  equations,  we  first  solve  for 
f^^^  from  Eq.  (7),  then  insert  f^^^  into  Eq.  (12)  to  calculate 
The  final  expression  for  the  dispersion  relation  is  then  obtained  by 
substituting  into  Eq.  (13)  and  carrying  out  the  r-  integration. 

Equation  (7)  can  be  solved  by  the  method  of  characteristics, 
namely,  by  integrating  it  along  the  unperturbed  (helical)  trajectories 
of  the  electrons.  Considerable  algebra  is  involved  in  this  integra- 
tion; however,  the  procedures  are  straightforward  and  standard.  Here 
we  simply  present  the  final  result. 


where  Q - eB  /me,  ip  and  6 are  the  polar  angles  of  p 

GO  ~ X 

is  the  Bessel  function  of  order  s,  y = [1  + (p^  + 
Gssi(x)  = J5+51 (x)dJj(x)/dx.  Inserting  Eq.  (14)  into 
substituting  the  resulting  equation  into  Eq.  (13),  we 


and  r,  respectively, 
and 

Eq,  (12)  and 
obtain 
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S ' -1  y.  W 


J-  , 2 2 

- "n 
c 


4Tre‘ 


mc^r  ^(a  r ) 

w 0 ' n ^ ^ , 
s , s - 


2 / rJj(^.pr)J^.('.^r)dr 


2tt 


.y*  y*  dp^y  d.j>  Jexp  [-i(s'-l)i|;]  -exp  [-i(s'  + l 


)'>]! 


ss 


' \ / 


Y'^-k  p /m-sQ„ 
z z e 


L . i ^ 

\ / p,  '<p. 


k 9f. 


Ym  '^p 


(15) 


where  'p  ^ ip  - 0 


To  carry  out  the  r-integration  in  Eq.  (15),  we  need  to  first  construct 
a distribution  function  f^  from  the  constants  of  the  motion.  The 
distribution  function  consistent  with  assumption  (i)  is 


fo  = K 6(rL^-2c  Pg/eB^-r^^)  6(p^-p^^)  6(p^-p^q)/P^  , (16) 


where  6(x)  is  the  Dirac  delta  function,  p^^,  p^^  are  the  perpendicular 


and  parallel  momenta  of  the  electrons,  r = p /m  is  the  electron  Larmor 

2 L i>  6 


radius,  p^  = ymrp  cos  p - eB^r  /2c,  and  K is  a normalization  constant 
to  be  determined  from  the  equation 


/' 


f^  2Tirdrd  p = N 


where  N is  the  number  of  electrons  per  unit  length.  After  some  algebra, 
Eq.  (16)  can  be  rewritten  as 


8 


where 


fg  = — 2 — + 6('1'-tt+'1>q)]  '5(p^-Pj^o)  '^(Pz'Pzo^ 

2tt  p 

i 


•S(r-rj)  S(r2-r) 


Jq  = sin"^  [(r^  + rj-  - / 2rr|^]  , 


(l  , X > 0 

" |0  , X < 0 

-■l  ■ I ■'o  ■ ■'l  I 


•"2  " '"o  ""l  * 


Substituting  Eq.  (17)  into  Eq.  (15),  we  can  carry  out  all  the  integra- 
tions exactly  and  thereby  obtain  the  dispersion  relation 


"7  - "z  - “n 
c 


-4v  n..^-k/c^) 


(».-l<,v,„)Q,(a„ro.»„rL) 

(»  - - sa^) 


2 222  1/2 

where  Y„  - [1  ♦ (P  * Pj„  )/»i  c ] ' , 6^^  • P^„/Y„i"c.  = Pzo'V' 
^c  " ^e^’*^o’  ~ NeVmc^  is  a dimensionless  beam  density  parameter,  and 
the  functions  and  Qg  are  defined  as 
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(19) 


= [J3  (x)J'(y)]2 


Q^(x,y)  H 2H^(x,y)  + yj;{y)J^(y) ^j;(x)(l  + sVx^)  + [j;(x)]-} 

+ 2s“J^{x)J^(x)J^(y)[yJ^{y)  - J3(y)]/xy  (20) 

B.  Normalization  and  Beam  Frame  Representation 

For  generality,  the  following  normalization  procedures  will  be 
introduced  through  which  the  waveguide  radius  r^  is  scaled  out  of  the 


dispersion  relation: 

(i)  Length  normalized  to  r^,  (21) 

w 

(ii)  Frequency  normalized  to  c/r  , (22) 

(iii)  Momentum  normalized  to  me,  and  (23) 

2 

(iv)  EM  fields  normalized  to  me  /er  . (24) 

w 


Furthermore,  we  will  present  the  results,  both  analytical  and  numerical, 
in  the  reference  frame  in  which  the  beam  axial  velocity  vanishes.  The 
advantage  for  introducing  the  normalization  procedures  and  the  beam 
reference  frame  is  that  the  numerical  data  to  be  presented  will  be 
applicable  to  any  waveguide  radius  and  any  beam  axial  velocity  through 
appropriate  transformations.  In  Appendix  A,  we  have  included  the 
transformation  formulae  for  conversion  from  beam  to  lab  frame  quantities 
and  from  normalized  to  physical  units. 

To  maintain  a clear  distinction  between  lab  frame,  beam  frame, 
normalized  and  physical  quantities,  the  following  notation  has  been 
adopted: 

(i)  Beam  frame  quantities  are  denoted  by  primed  symbols 
(e.g.,  y')»  laboratory  frame  ouantities  are  unprimed 
(e.g..  Y°). 

(ii)  Normalized  diiriens ionless  quantities  are  denoted  by  a bar 
(e.g.,  w).  Unbarred  symbols  represent  either  physical 
Quantities  (e.g.,  w)  or  naturally  dimensionless  ones 
(e.g.,  Yq)- 
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In  beam  frame,  the  normalized  form  of  the  dispersion  relation 
[Eq.  (18)]  can  be  written 


-■2 

(ts 


y'J^(x  ) 
0 0 n' 


(w'  - sYi'^)^ 


(25) 


Note  that  the  primes  in  Eqs.  (19)  and  (20)  represent,  respectively, 

the  first  and  second  derivatives  of  the  Bessel  function  with  respect 

to  its  argument.  In  Eq.  (25)  and  elsewhere  in  this  paper,  the  primes 

indicate  beam  frame  quantity.  In  Eq.  (25),  and  are  frame 

independent  quantities.  It  can  be  easily  shown  that  the  normalized 

cutoff  frequency  of  the  TE^^  mode,  is  also  frame  independent  and 

numerically  equal  to  x . 

n 


C . Threshold  Condition  for  Instability 

On  the  right  hand  side  of  Eq.  (25),  the  first  term  is  the 
source  term  of  the  instability.  However,  instability  occurs  only  when 
exceeds  a certain  threshold  value.  To  find  this  threshold  value, 
we  expand  Eq.  (25)  about  the  intersecting  point  of  the  waveguide 

characteristic  curve 

—'2  — ? 

■ ''z  ■ n " ° 

and  the  beam  characteristic  curve 


0) 


(27) 


n 


where 

U)'  = + 

0 


is  chosen  so  that  the  two  curves  do  intersect. 
" ^zo’  ^c  = ^co 


Substituting 


into  Eq.  (25),  we  obtain 


2v's'r2'  0 2v'x^H 

.-  3 CO  s , n s iO 

Au)  - o Aoj  + 


Y u)'J^(x  ) 
0 0 o'  n' 


Y^w' J^(x  ) 
'0  0 n 


= 0 


(28) 


The  condition  for  Eq.  (28)  to  have  complex  solutions  is 


3 > 


8v' 


1/4  (sn;  0 )^/'‘ 


co^s' 


n s 


1/2 


(29) 


Equation  (29)  can  be  solved  for  the  threshold  velocity  for  unstable 
interactions.  Note  that  and  are  both  functions  of  3'  • 

5 S xO 

D.  Optimization  of  the  Beam  Guiding  Center  Position  r^ 

The  function  in  the  source  term  of  the  instability  indicates 
the  strength  of  beam-wave  coupling  with  respect  to  r^.  A physical 
interpretation  of  can  be  found  in  Ref.  15.  Figure  3 plots 
^s^*n^o’^n^L^  ^ function  of  its  arguments.  We  note  that  is  a 

nonmonotonic  function  of  r^  with  an  infinite  number  of  peaks.  The 
first  peak  of  has  the  highest  amplitude,  thus  the  value  of  which 
corresponds  to  the  first  peak  of  is  the  optimum  choice  for  F^. 

Table  I lists  the  optimum  values  of  F^  for  various  combinations  of  s 
and  n.  These  values  will  be  used  later  in  selecting  parameters  for 
number! cal  computations. 


E.  Optimization  of  the  Wave  Number 

For  gyro-TWA's  operating  at  the  nonfundamental  cyclotron  harmonics 
the  beam  energy  (usually  tens  of  keV)  is  well  above  the  threshold 
given  by  Eq.  (23),  hence  the  second  term  of  Eq.  (25)  can  be  neglected. 
Again  we  expand  Eq.  (25)  about  the  point  However,  we  will 


now  expand  in  terms  of  three  variables  - w' , k^.  and 
just  one.  Substituting 


= u) ' + Aw' 


instead  of 


4o  ^ 


into  Eq.  (25)  and  neglecting  the  second  term  on  the  right  hand  side, 
we  obtain 


— * ^zn  — * — ' — * 2 

Aw  - — T"  Ak  (Aw  - SASl^) 
— z c 

w 


I p 

'o  0 n 0 


In  Eq.  (30),  Aw  = Aw  + iAw.  is  a complex  number,  while  all  other 

' ' I 

quantities  are  real.  Differentiating  Eq.  (30)  with  respect  to  Ak^, 
we  obtain 


d(Aw  ) 
d(Ak^) 


- SAfij.) 


Wq(3Aw  - sAfi^  - 2ir^QAl^2^‘^0^ 


I I 

Equation  (31)  shows  that  if  k„  = 0,  the  growth  rate  Aw.  is  a maximum 

I I 

regardless  of  the  values  of  Alc^  and  Afi^,  i.e.,  Aw^  remains  a maximum 


- 
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<»v.; 


▼ 


1 


/w!  = /3 


'2 

Vo'-n) 


1/3 


(35) 


* 


I 

» 

I 

f • 
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Equation  (35)  is  an  analytical  expression  of  the  maximum  growth  for 
the  usual  case,  where  the  beam  energy  far  exceeds  the  threshold 
condition. 


Inserting  Eq.  (23)  into  Eq.  (29),  we  obtain  the  threshold 
condition  for  the  case  of  grazing  intersection 


8v' 


27y„J^  (x„) 
0 0 n 


1/4  Q, 


3/4 


, 1/2.  1/2 

^n  ^s 


(36) 


G . Numerical  Solutions  of  the  Dispersion  Relation 

I 

Figure  5 shows  the  growth  rate  ‘ as  a function  of 
for  the  first  four  cyclotron  harmonics.  The  beam  parameters  used 
were  v'  = 0.002  and  = 1.1,  which  characterize  the  electron  gun 
designed  for  the  NRL  gyrotron  program.  In  Fig.  4,  the  magnetic  field 
was  chosen  so  that  X '=  1 (i.e.,  the  grazing  intersection  case).  Also 
the  radial  eigenmode  number  n was  set  equal  to  the  cyclotron  harmonic 
number  s so  that  the  optimized  r^  as  found  from  Table  I (the  numbers 
marked  by  "*")  indicates  a beam  position  centrally  located  between 
the  waveguide  axis  and  wall.  Henceforth,  we  will  always  let  n equal 
s in  the  numberical  computations  and  the  corresponding  r^  in  Table  I 
will  be  used.  However,  the  analytical  scaling  relations  [e.g.,  Eqs.  (34), 
(35),  (45),  and  (50)]  can  be  used  to  scale  the  numerical  results  to 
other  combinations  of  n and  s.  The  growth  rate  plots  in  Fig.  4 are 

I I 

typical  solutions  of  Eq.  (25).  One  notes  that  peaks  at  IT  = 0 
and  decreases  as  either  |k^|  or  s Increases.  Information  concerning 
power  gain  and  bandwidth  can  also  be  deduced  from  these  plots  (discussed 
next) . 
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H. 


Power  Gain  and  Bandwidth 


Let  G be  the  total  power  gain  of  a gyro-TWA  system  which  is 
characterized  by  a linear  growth  rate  and  interaction  length  L. 

At  the  linear  growth  rate,  the  wave  power  amplification  is  given  by 

P + = P.  exp  (2(jj.t)  (37) 

out  in  ^ ' 1 ' 

where  t h L/v^^,  v^^  is  the  beam  axial  velocity  and  P.^,  P^^^^  are, 

respectively,  the  input  and  output  wave  powers.  In  terms  of  G,  P^^ 

and  P ^ are  related  through 
out 

P . = . (38) 

out  in 

From  Eqs.  (37)  and  (38),  we  obtain  the  total  power  gain 
G = 20  (jj^x/JlnlO 

= 8.7  (j.T  dB  (39.1) 

and  the  power  gain  per  unit  length 

g = 8.7(o^/V2Q  (39.2) 

To  find  the  bandwidth,  we  let  u)^  and  J]  be  the  growth  rates 
which  yield  the  maximum  power  and  half  of  the  maximum  power,  respectively, 
and  let  g"*  and  G^  be  the  corresponding  power  gains.  By  definition  g"" 
and  G^  are  related  through 


G^  = g"*  - 10  log  2.  (40) 


R = (oj/o)'!'  = G^g'"=  1 - 10  log  2/g'". 


(40) 


This  ratio  determines  the  half  power  points  in  the  plot  (see 

Fig.  4)  and  therefore  the  width  aIT^  between  the  half  power  points. 
Through  a Lorentz  transformation,  we  obtain  the  lab  frame  bandwidth 

A6j/6j  ts 

where  v^^/c. 

As  an  example,  we  assume  G = 20,  hence  Eq.  (40)  gives  R = 0.85. 
From  Fig.  4 and  Eq.  (41),  we  find  that  the  bandwidth  for  s = 1,  2,  3, 
and  4 are,  respectively 

Aoj/oj  = 11.19%,  7.44%,  5.11%,  and  3.91%. 


I . Sensitivity  to  X' 

We  consider  the  special  case  of  grazing  or  near  grazing  inter- 
sections (X'  < 1,  = 0).  Thus,  substituting  = 0 and 

into  Eq.  (30),  we  obtain 

Aw'(A(:r'  - sAii--^)^  = - 2v'x^H^3j^g^/YgJo^(xJ.  (42) 


From  Eq.  (42),  we  obtain  a condition  for  unstable  interactions, 

2 


sAil^  > -3 


v'Vs^lo 


(x„) 

0 0 n' 


1/3 


(43) 


Since  X'  = s(n'  „ + An'  )/x„  - 1 + sAfi  Vx^,  Eq.  (43)  can  be  written  as 
CO  c n c n 


X'  > X' 


min 


(44) 
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where  X’  . is  defined  as 
mi  n 


’^'min  =1-3 


v'H  P' 

S 10 


^o’^nV(^n) 


1/3 


(45) 


In  Sec.  IV,  we  show  that  by  lowering  X'  from  unity  toward 

the  efficiency  increases  sharply.  However,  it  can  be  shown  from 

Eq.  (39)  that  X^^.^  is  generally  very  close  to  unit,  which  implies  that 

efficiency  is  a very  sentitive  function  of  X'.  Thus  in  order  to 

optimize  the  efficiency  by  adjusting  X',  one  should  be  able  to  accurately 

control  X'  between  1 and  X’  . . This  point  will  be  discussed  quanti- 

mi  n 

tatively  in  Sec.  V for  specific  cases. 
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APPROXIMATE  ANALYTICAL  SCALING  RELATION  FOR  THE  EFFICIENCY 


III . 

From  the  foregoing  analysis,  it  is  clear  that  a gyro-TWA  can 
ooorate  on  various  combination  of  s and  n and  involve  numerous  parameters. 
By  scaling  out  the  waveguide  radius  and  working  in  the  beam  frame,  we 
have  generalized  the  results  considerably.  Still  the  remaining  parameters 
are  too  many  for  numerical  optimization.  This  is  especially  true 
for  cases  involving  nonlinear  calculations.  In  this  regard,  general 
scaling  relations  prove  invaluable.  In  the  linear  analysis,  we  have 
already  shown  now  the  growth  rate  [Eg.  (35)],  instability  threshold 
[Eg.  (36)],  and  sensitivity  to  X’  [Eg.  (45)]  scale  with  parameters 
such  dS'J,y^\  n,  and  s.  Thus,  from  the  numerical  data  on  one  set 
of  parameters,  one  can  deduce  fairly  accurate  information  for  other 
parameters.  Similarly,  for  the  nonlinear  analysis  that  follows,  we 
will  complement  the  numerical  approach  with  analytical  scaling 
relations  so  that  discrete  numerical  data  points  can  be  interpolated 
and  extrapolated  to  cover  essentially  the  entire  parameter  space. 

In  Sec.  II,  we  have  shown  that  the  following  condition  holds 
at  the  onset  of  the  instability, 

^'r  ■ ^^e^’*''o  ""  ^“'r’ 

where  Aoj  is  given  by  Eg.  (28)  for  the  case  of  grazing  intersection 

— ^ ft 

and  = 0.  Recent  nonlinear  analysis  has  shown  that  a similar 

condition  holds  at  the  saturation  state, 

where  <y’^>  is  the  average  y'  of  all  electrons  at  saturation.  Eg.  (46) 
minus  Eg.  (47)  gives 


si^e  ■ ^0  2A(o*  . (48) 


A 


From  Eq.  (48)  and  assuming  A6d'^  <<  x^,  we  obtain  an  approximate 
estimate  of  the  efficiency  n'  (in  beam  frame). 


- 1 'W^c 

0 

For  the  case  of  grazing  intersection  (X'=l, 
(49)  give 


(49) 

Xn  = sf2^),  Eq.  (34)  and 


n'  = 


2v'  H.  8'/ 
s iO 


y' 

'0  n o n 


1/3 


(50) 


Equation  (50)  shows  how  n ' scales  with  y^,  s,  n,  and  v'.  The 
scaling  with  v'  and  n is  particularly  simple,  namely. 


n ’ ct  V 


1/3 


(51) 


and,  for  the  optimum  beam  positions  given  in  Table  I, 

n-a  [x^Jo(x^)]'2/^  (52) 

For  n = 1,  2,  3,  4,  and  5,  the  right  hand  side  of  Eq.  (52)  equals 
0.75,  0.61,  0.54,  0.49,  and  0.45,  respectively.  Thus,  using  lower 
radial  eigenmode  numbers  generally  results  in  higher  efficiency. 
However,  we  note  that  for  the  same  output  wave  frequency,  a lower 
radial  eigenmode  number  implies  a smaller  waveguide  radius. 
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IV. 


SINGLE  WAVE  SIMULATION 


A single  wave  numerical  simulation  code  has  been  developed  to 

study  the  saturation  mechanisms  and  calculate  the  energy  conversion 

efficiency.  This  code  is  a cylindrical  version  of  the  Cartesian 

7 8 

code  originally  developed  by  Sprangle,  Drobot,  and  Manheimer  ’ . 

For  a detailed  description,  see  Ref.  8.  The  simulation  shows  two 

g 

different  saturation  mechanisms  : depletion  of  free  energy  and  loss 
of  phase  synchronism.  The  first  mechanism  dominates  when  the  beam 
energy  is  only  slightly  above  threshold.  Saturation  occurs  as  soon 
as  the  beam  loses  a small  amount  of  energy  and  the  system  becomes 
linearly  stable.  This  mechanism  can  lead  to  saturation  even  in  the 
linear  stage.  The  second  mechanism  dominates  when  the  beam  energy 
is  well  above  the  threshold.  Saturation  occurs  because  an  average 
electron  loses  so  much  energy  that  its  relativistic  cyclotron 
frequency  $7^  no  longer  matches  to  favor  unstable  interactions. 

As  found  in  Ref.  8,  both  mechanisms  are  important  for  the  fundamental 
cyclotron  harmonic  because  the  beam  threshold  energy  is  not  negligible. 
However,  for  nonfundamental  cyclotron  harmonics,  we  found  that  the 
threshold  energy  is  usually  so  low  that  the  first  mechanism  can  be 
disregarded.  This  is  a basic  difference  between  the  fundamental 
and  nonfundamental  cyclotron  harmonic  interactions  as  far  as  satura- 
tion is  concerned.  The  scaling  relations  derived  in  Sec.  II  and 
Sec.  Ill  by  neglecting  the  threshold  condition  would  thus  be  generally 
valid  for  nonfundamental  cyclotron  harmonics.  In  fact,  all  the 
simulation  data  which  follows  agrees  with  the  scaling  relation  [Eg. 
(50)]  to  within  10%.  In  other  words,  from  one  single  data  point,  all 
other  data  points  can  be  predicted  with  an  accuracy  better  than  90%. 
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As  a cross  check  of  the  code  the  linear  frequency  and  growth 
rate  obtained  from  Eq.  (25)  have  been  compared  to  the  linear  stages 
of  the  simulation.  The  agreement  was  very  good.  Also,  conservation 
of  energy  was  monitored  in  the  code  and  never  deviated  by  more  than 
0.1%. 

There  are  seven  free  parameters  to  be  specified  in  making  a 
simulation  run,  i.e.,  s,  n,  k^,  v',  and  X'.  For  reasons 
described  in  Sec.  II,  we  will  let  s and  n be  equal.  Linear  analysis 
of  Sec.  I has  shown  that  maximum  beam-wave  coupling  occurs  when 
and  r^  assumes  the  values  in  Table  I.  Thus  there  are  four  remaining 
parameters:  s,  u',  Yq>  and  X'.  We  will  restrict  our  study  to  the 
third  and  fourth  cyclotron  harmonics  and  let  u'  = 0.002  and  0.005. 

The  beam  energy  Yq  will  be  varied  from  Yq  1 .04  to  1 .5  (i .e. , 20  keV 
to  250  keV  in  the  beam  frame).  Finally,  the  specification  of  X'  is 
based  on  the  following  consideration.  Figure  6 shows  two  typical 
plots  of  the  efficiency  (n ' ) and  linear  growth  rate  (cop  vs  X'  for 
v'  = 0.005,  Yq  = 1-1.  and  s = 3,  4.  It  can  be  seen  that  both  n'  and 
ooj  are  sensitive  functions  of  X'.  However,  as  shown  in  Fig.  6,  when 
n'  increases,  co^  always  decreases.  High  efficiency  achieved  at  too 
low  a growth  rate  is  clearly  undesirable.  Thus  some  compromise  must 
be  made  between  the  growth  rate  and  the  efficiency  in  our  selection 
of  X'.  We  have  decided  to  select  two  values  of  X':  one  at  the  maximum 
growth  rate  (X'=l),  the  other  at  half  of  the  maximum  growth  rate 
[X'=X|:^,  where  X^  is  to  be  found  from  Eq.  (25)].  The  efficiency  calcu- 
lated at  X'=X,'  will  be  referred  to  as  the  optimized  efficiency, 
h 

although  it  could  still  be  increased  if  one  can  tolerate  a lower  growth 
rate.  Figure  7 plots  n'  versus  Yq  s=3  and  the  selected  values 
of  v’  and  X'.  One  notes  that  n'  is  a monotonical ly  decreasing  function 
of  Yq-  We  did  not  continue  the  curves  through  the  low  y^  region, 
because  the  growth  rates  become  too  low  in  that  region  and  simulation 
runs  become  exceedingly  expensive.  The  points  at  Yq=1-06  and  1.1 
marked  by  1,  2,  3,  and  4 have  been  selected  for  experimental  designs 


(see  Sec.  V)  because  they  are  close  to  the  peak  efficiency.  Lower 
beam  energies  will  result  in  even  higher  efficiencies,  but  beam  power 
and  growth  rate  will  both  be  lower.  Figure  8 plots  the  same  curves 
for  s=4.  In  contrast  to  the  case  of  s=3,  here  the  efficiency  is 
generally  lower  and  has  a peak  at  y 1.14.  Such  differences  are 
also  predicted  analytically  [Eq.  (50)].  The  points  at  = 1-10 
and  1.14  marked  by  5,  6,  7,  and  8 have  been  chosen  for  experimental 
designs. 


V. 


EXPERIMENTAL  DESIGNS 


We  now  convert  each  of  the  selected  optimuni  or  near  optimum 
points  in  Figs.  7 and  8 into  an  experimental  design.  Parameters 
corresponding  to  these  points  are  beam  frame  normalized  quantities. 

To  convert  these  parameters  into  actual  design  parameters,  two  pieces 
of  information  are  required. 

(i)  the  beam  axial  velocity  (in  lab  frame),  which  can 
be  calculated  by  specifying  the  ratio  °‘t=''^o'^'^zo^ 
in  the  lab  frame.  We  let  a=1.5. 

(ii)  the  waveguide  radius,  which  is  determined  from  the 
desired  wave  frequency,  f=94  GHz,  [see  Appendix  A, 

Eq.  (A-4)]. 

Using  the  conversion  formulae  in  Appendix  A,  we  obtain  eight  p>^eliminary 
designs  for  a 94  GHz  Gyro-TWA  employing  the  third  or  fourth  cyclotron 
harmonic.  Tables  II  and  III  list  the  design  parameters.  Within  the 
context  of  the  present  model  and  assumptions,  the  main  difficulty  in 
achieving  the  projected  efficiency  appears  to  be  the  high  sensitivity 
associated  with  the  parameter  X'  (“B^/Yq).  As  shown  in  Fig.  6,  this 
parameter  has  to  be  tunable  with  an  accuracy  of  '^0.2%  for  s=3  and 
'^0.1%  for  s=4.  In  the  next  section  we  discuss  some  difficulties  not 
addressed  in  our  model  but  which  could  be  of  overwhelming  importance. 


VI. 


DISCUSSION 


The  designs  shown  in  Tables  II  and  III  all  represent  record 
microwave  powers  in  the  millimeter  frequency  range.  Despite  their 
attractiveness,  it  should  be  emphasized  that  these  designs  are  of  a 
preliminary  nature  because  questions  concerning  the  realization  of 
our  assumed  cold  beam  single  wave  model  remain  unanswered.  If  such 
an  idealized  model  is  not  realizable,  one  must  determine  how  thermal 
effect  and  competitive  mode  effect,  etc.,  alter  the  predicted 
efficiency.  Here  we  briefly  address  these  questions  in  a qualitative 
manner. 


where  is  the  axial  velocity  spread  and  Ay  is  the  energy  spread. 

On  the  basis  of  recent  simulation  studies  of  the  magnetron-type 
electron  gun^^,  we  have  found  that  condition  (53)  can  be  marginally 
satisfied  for  the  first  cyclotron  harmonic.  Since  condition  (53) 
becomes  more  restrictive  for  higher  harmonics,  it  appears  that 
further  improvement  of  the  electron  gun  design  is  required  to  generate 
a suitable  electron  beam.  Also,  thermal  effects  should  be  incorporated 
in  future  models  because  of  their  importance  in  higher  harmonic 
operations . 

We  have  examined  the  problem  of  competitive  modes  with  a 
linear  analysis^^.  All  intersections  of  the  waveguide  mode  charac- 
teristic curves  with  beam  cyclotron  harmonic  curves  were  found  to  be 
unstable  to  some  degree.  The  ones  travelling  backward  could  be 
classified  as  absolute  instabilities  and  might  present  the  greatest 


danger.  It  should  be  stressed  that  it  is  not  yet  clear  how  seriously 
these  inodes  will  affect  the  device  operation  and  further  analysis  is 
recommended.  On  the  other  hand,  this  problem  serves  to  underline 
the  advantage  of  beam  position  optimization  discussed  in  Sec.  II 
(item  D).  It  has  been  shown^^  that  through  such  optimization  the 
operating  mode  can  at  least  have  the  highest  growth  rate  as  compared 
with  all  the  competitive  modes. 

Finally,  we  point  out  that  the  beam  space  charge  effects  have 
not  been  examined  in  any  detail  in  the  present  report.  Presumably  it 
will  reduce  the  efficiency  by  introducing  beam  thermal  and  energy 
spread. 
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APPENDIX  A 


CONVERSION  FORMULAE 


For  generality,  numerical  data  are  presented  for  beam  frame 
normalized  quantities.  In  this  appendix,  we  tabulate  the  formulae 
needed  to  convert  beam  frame  normalized  quantities  to  lab  frame  physical 
quantities.  Column  2 of  Table  IV  converts  the  beam  frame  normalized 
quantities  in  column  1 to  lab  frame  normalized  quantities.  The  infor- 
mation needed  for  the  conversion  is  the  axial  velocity  v^^  (in  lab 
frame)  which  defines  the  quantity  in  column  2, 


ir,  = II 


Formulae  presented  in  column  2 are  based  on  the  Lorentz  transform. 
Derivations  are  obvious  except  for  the  following  items. 


is 


(i)  n = -1)  (A.l) 

To  derive  Eq.  (A.l),  we  note  that  the  wave  energy  per  unit  length 

"2  2 
kbV 

Wf  = ^ 

c 


where  K is  a geometrical  factor,  is  the  axial  wave  magnetic  field 
amplitude,  and  co^  is  the  waveguide  cutoff  frequency. 

The  beam  kinetic  energy  per  unit  length  is 
^b  ” 

where  N is  defined  in  the  text. 


The  lab  frame  efficiency  (n)  can  therefore  be  written 

^2  2 2 
W,  KB^mc  w 

^ I ■ Z I— ■. — ■ 

‘''b  8TTa,2^(Yo  -1) 


(A. 2) 
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In  Eq.  (A.2),K,§^  and  u)^  are  all  frame  independent  quantities. 


► 

s 


f 
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The  beam  frame  efficiency  (n')  is  also  expressed  by  Eq.  (A. 2) 
with  0),  N,  and  y replaced  by  w',  N',  and  y'  respectively.  Thus 


ri'  (y  -1) 


(A.3) 


Since  w ' ^ substitution  of  these 

relations  into  Eq.  (A.3)  yields  Eq.  (A.l). 


(ii)  “"'i  ^ 


In  the  beam  frame,  k'  = 0.  Hence  all  quantities  vary  as  exp 

(-iw't').  To  transform  this  expression  into  the  lab  frame,  we  use 

2 

the  relation  t'  = y^{t  “V^qZ/c  ) and  obtain 

2 

exp  ( -iu't')  = exp  I -i(u)„'  + iw")  Y,(t  -v  z/c  )] 


= exp  1 -iY^w^'t  + i- 


Y_w  'V. 

' ^ i/» 


Z + - 


Y,a)/v. 


z . 


Thus,  the  exp  (-iw^t')  dependence,  viewed  from  the  lab  frame,  has 
the  following  wave  frequency  (w^)  and  wave  number  (k^,). 


0)  = Y (A) 

r r 


''z  = 


and,  if  one  follows  a beam  segment  (i.e.  let  z = v^^t),  one  observes  a 
growth  rate  given  by 


-♦••Ff 


Column  3 converts  the  lab  frame  normalized  quantities  in  column  2 
to  lab  frame  physical  quantities.  The  conversion  formulae  are  based  on 
Eqs.  (21)-(23)  in  the  main  text.  The  information  needed  at  this  step  is 
the  waveguide  radius  r , which  can  be  determined  from  the  desired  wave 

W 

frequency  through 


i^w  = u)^c/2TTf,  (A.4) 

where  is  the  calculated  wave  frequency  (normalized)  and  f is  the  de- 
sired wave  frequency  in  Hz. 


From  the  lab  frame  physical  quantities  in  column  3,  we  can  derive 
the  beam  current  beam  power  (Pjj)t  and  the  maximum  output  wave 


power  (P^). 

U O' 

Ijj  = vmc^v^^/e  = 1.707  x 10^ 

(A.5) 

^b  = 

Ib(Yo  = 8.535  x 10^(Yo 

(A.6) 

^w  = ^Pb- 

(A.7) 

I 


s. 
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Table  I.  Values  of  7^  for  which  absolute  maximum 

with  respect  to  the  argument  ^^7^-  The  numbers  marked  by  * have  been  used 
to  obtain  the  data  presented  in  Fig.  5 and  also  in  simulation  runs. 


Table  II.  Efficiency  optimized  design  parameters  for  a 94  GH^  gyro-TWA 
opearting  at  the  third  cyclotron  harmonic.  Data  no  refers  to  the  points 
marked  on  Fig.  7. 


data 

no 

1 

2 

3 

4 

V 

(XIO^) 

2.05 

2.08 

5.12 

5.19 

Wb 

(keV) 

42.80 

70.82 

42.80 

70.82 

'b 

(Amp) 

7.55 

9.48 

18.87 

23.70 

n 

(%) 

12.73 

11.57 

15.41 

14.93 

Pb 

(kW) 

323.1 

671.5 

807.8 

1679.0 

(kW) 

41.2 

77.7 

124.5 

250.7 

»o 

(kG) 

11.55 

11.81 

11.54 

11.78 

'z 

(cm"^) 

4.25 

5.27 

4.25 

5.27 

*^w 

(mm) 

5.285 

5.356 

5.285 

5.351 

*"0 

(mm) 

2.182 

2.212 

2.182 

2.209 

'■l 

(mm) 

0.520 

0.663 

0.520 

0.664 

"1 

(mm) 

1.663 

1.549 

1.662 

1 .545 

(mm) 

2.702 

2.874 

2.702 

2.873 

Vjjo/C 

0.324 

0.401 

0.324 

0.401 

0.216 

0.268 

0.216 

0.268 

9 

(dB/cm) 

1.167 

1.583 

1.844 

1.855 

Table  III.  Efficiency  optimized  design  parameters  for  a 94  GH^  gyro-TWA 
operating  at  the  fourth  cyclotron  harmonic.  Data  no  refers  to  the  points 
marked  on  Fig.  8. 


data 

no 

5 

6 

7 

8 

V 

(XIO^) 

2.08 

2.10 

5.19 

5.25 

(keV) 

70.82 

98.49 

70.82 

98.49 

^b 

(Amp) 

9.48 

10.93 

23.70 

27.32 

n 

(%) 

6.67 

6.96 

8.87 

9.25 

Pb 

(kW) 

671.5 

1076.0 

1679.0 

2691.0 

Pw 

(kW) 

44.8 

74.9 

149.0 

248.8 

«o 

(kG) 

8.878 

9.085 

8.872 

9.070 

^z 

(cm'b 

5.27 

6.00 

5.27 

6.00 

(ran) 

7.017 

7.093 

7.011 

7.093 

"o 

(ran) 

2.800 

2.831 

2.798 

2.831 

(ran) 

0.881 

1.028 

0.882 

1.030 

"l 

(ran) 

1.919 

1.802 

1.916 

1.801 

•^2 

(ran) 

3.681 

3.859 

3.680 

3.861 

V^/C 

0.401 

0.457 

0.401 

0.457 

V^/C 

0.268 

0.305 

0.268 

0.305 

g 

(dB/cm) 

0.738 

0.779 

0.930 

1.130 

I 


Table  IV.  Formulae  for  transformation  from  beam  frame  into  lab  frame 


and  for  converting  normalized  quantities  into  physical  quantities. 


1 

Beam  Frame 

Lab  Frame 

Lab  Frame 

Normalized  Quantities 

Normalized  Quantities 

Physical  Quantities 

j 

1 

V 

v=Y2v' 

v(dimensionless) 

Y' 

Y=Y2Y' 

Y(dimensionless) 

n' 

n=n'Y2(Y'  -1)/(Y2Y'  -1) 

n(dimensionless ) 

4 (=0) 

'‘z  ■ Vzc,"r 

“r  = ^z^r 

'^r  = 

'V 

W.  = w1/y_ 

1 1 z 

o^i  = w.c/r^ 

a' 

n = 

Q = n'c/r, . 

e 

e e 

e e w 

\ 

•"l  " '"l’"w 

r' 

r = r' 

Y = r„r 

0 

0 0 

0 0 w 

I 

f 

I 

I 

♦ 

1 

t 

I 
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Instantaneous  positions  of  two  electrons  circulating  in 
the  same  cyclotron  orbit,  but  separated  by  180°  in  phase 
angle.  E is  the  wave  electric  field  and  (pointing 
toward  the  reader)  is  the  applied  magnetic  field. 


Cross-sectional  view  of  the  electron  cyclotron  maser 
model.  The  applied  magnetic  field  (not  shown)  points 
toward  the  reader.  The  electrons  are  monoenergetic  and 
all  have  the  same  Larmor  radius  r^^.  Guiding  centers  of 
all  electrons  are  uniformly  distributed  on  the  circle  of 
contant  radius  r.. 
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Fig.  4. 


The  vacuum  wave  guide  characteristic  curve  and  the  beam 
characteristic  curve  in  the  plane  (plotted  in  the 

beam  frame). 
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Fig.  5. 

Growth  rates  (wp  versus  for  the  first  four  cyclotron 
harmonics  as  computed  from  Eq.  (25).  Parameters  used  for 
this  figure  were  v'=0.002,  =1-1  (50  keV),  X'=l,  and 

Yp=0.48,  0.43,  0.41,  0.40  for  s=l , 2,  3,  4 respectively. 
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